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ABSTRACT 

The production factor, or broad band averaged cross-section, for solar wind charge-exchange with 
hydrogen producing emission in the ROSAT | keV (R12) band is (3.8 db 0.2) x 10 -20 count degree -2 
cm 4 . This value is derived from a comparison of the Long-Term (background) Enhancements in the 
ROSAT All-Sky Survey with magnetohysdrodynamic simulations of the magnetosheath. This value is 
1.8 to 4.5 times higher than values derived from limited atomic data, suggesting that those values may 
be missing a large number of faint lines. This production factor is important for deriving the exact 
amount of | keV band flux that is due to the Local Hot Bubble, for planning future observations 
in the \ keV band, and for evaluating proposals for remote sensing of the magnetosheath. The 
same method cannot be applied to the | keV band as that band, being composed primarily of the 
oxygen lines, is far more sensitive to the detailed abundances and ionization balance in the solar wind. 

We also show, incidentally, that recent efforts to correlate XMM-Newton observing geometry with 
magnetosheath solar wind charge-exchange emission in the oxygen lines have been, quite literally, 
misguided. Simulations of the inner heliosphere show that broader efforts to correlate heliospheric 
solar wind charge-exchange with local solar wind parameters are unlikely to produce useful results. 
Subject headings: X-rays: diffuse emission 


1. INTRODUCTION 
1.1. Motivation 


Solar wind charge-exchange (SWCX) occurs when a 
high charge state ion in the solar wind charge-exchanges 
with a neutral; the resultant ion is in an excited state and 
in the transition to the ground state a photon is produced 
in either the soft X-ray or extreme ultraviolet. Since the 
solar wind is temporally variable in density, speed, ele¬ 
mental abundance, and ionization fraction, the observed 
SWCX emission varies as a function of time, look di¬ 
rection, and observatory location. The charge-exchange 
spectrum has the same transitions, from the same ions, 
with nearly the same strengths, as a purely recombining 
plasma, and thus makes a strong contribution to the lines 
used for astrophysical plasma diagnostics ( See, for exam¬ 
ple, the discussion in Wargelin et al.|2008). Because the 
SWCX emission will be smooth over the EOV of any re¬ 
cent, current, or near future X-ray instrument, and since 
it has no unique spectral signatures, it is a particularly 
problematic foreground to remove from observations of 
any diffuse emission that also fills the FOV, such as that 
expected from the Local Hot Bubble, the Galactic bulge 
and halo, and the Warm-Hot Intergalactic Medium. It 
is also problematic for many observations of clusters of 
galaxies, nearby galaxies, local supernova remnants, and 
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super bubbles. 

Indeed, the uncertainty in the strength of the SWCX 
emission in the ROSAT | keV band has led to significant 
controversy over the existence of the Local Hot Bubble 
(LHB). The clearest evidence for the soft X-ray emission 
from the LHB is significant emiss ion in front of nearby 
(^ 60 pc) molecular clouds (e.g., Snowden et al. 1993). 
The distribution of hot gas was derived from an all- 
sky anti-correlation analysis fSnowden et al.|1998) which 
showed a rather hourglass shaped distribution where the 
waist is in the Galactic plane, and the bulges extend to¬ 
wards the Galactic poles. If, however, the bulk of the 
emission attributed to the LHB were actually due to 
SWCX, then the LHB might not exist, and the distri¬ 
bution of the re mainin g hot gas could be explained in 
other ways (|Welsh et al.|1999). Various other data, such 
as the pressure differenc e between the hot gas and the 
local interstellar clouds penkh^[2009), the distribution 
of O VI0 emission from t he LHB-cool cloud interfaces 
(Welsh & Lallement 2QQ8| among others), the shape of 
the local cavity (|Sfeir et al. 19 99]) , and the magneti c field 


in the wall of the LHB (Andersson & Potter 


2006) were 


variously deployed to argue both for and against the ex¬ 
istence of the LHB. While it was generally recognized 
that some part of the soft X-ray emission previously at¬ 
tributed to the LHB is indeed due to SWCX, estimates of 
the fraction of the emi ssion due to SWCX varied greatly 
(compare, for ex ample Koutroumpa et ah]|2009 ; Robert¬ 
son et al. [2009). Further contusion arose because the 
more reliable calculations for the SWCX emission were 
made in the | keV band, which is dominated by a few 
strong lines, while the bulk of the LHB emission is in 
the \ keV band where SWCX calculations are both diffi- 


6 We will use O VII to refer to the emission line, while we will 
use 0+ 6 to refer to the parent ion. 
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cult and highly uncertain. That uncertainty arises from 
the lack of reliable charge-exchange cross-sections for the 
multitude of faint lines that form the bulk of the | keV 
band emission. 

In the intervening years, observation al ver ific atio n of 
calculations of SWCX emission (e.g., Koutroumpa 2012) 
has focussed on the O VII and O VIII lines because 
those lines are accessible to Chandra , XMM-Newton , 
and Suzaku , while the j keV band is not. While such 
work has explored issues concerning the distribution of 
the neutral material with which the solar wind charge- 
exchanges and the structure of the solar wind itself, 
progress on the validation of our calculation of the SWCX 
emission in the | keV band has been minimal. While the 
number of measured or calculated charge-exc hange cross- 
sections has grown (|Betancourt-Martinez et al.|2014j and 
several other studies in preparation), these do not yet 
form a sufficient database with which to calculate reli¬ 
ably the SWCX emission in the | keV band. However, 
even if all of the contributing charge-exchange cross sec¬ 
tions and branching ratios were well known, modeling of 
the SWCX emission would still be problematic because of 
the lack of sufficiently detailed abundance and ionization 
state information for the solar wind. 


The DXL sounding rocket observation (Galeazzi et al. 
2014) of the He focussing cone (the region where the 
Sun gravitationally focuses the interstellar He flowing 
through the solar system) in the Wisconsin C-band al¬ 
lows a more empirical approach; the ratio of the observed 
X-ray emission to the model emission measure provides 
a broad-band abundance-weighted cross-section, or pro¬ 
duction factor. Although the Wisconsin C-band is not 
exactly the ROSAT | keV band, they are sufficiently 
similar that scaling from one to the other does not intro¬ 
duce significant uncertainties. However, DXL only pro¬ 
vides the production factor for SWCX with He (the IS 
H being strongly depleted near the He focussing cone), 
whereas a comparable emission is expected from SWCX 
with H, both from the interstellar H flowing through the 
solar system and the exospheric H in the Earth’s mag¬ 
netosheath. Here we derive a ROSAT \ keV band pro¬ 
duction factor for SWCX with H that requires no more 
than existing data and the application of relatively off- 
the-shelf magnetohydrodynamic (MHD) simulations. 

This work has implications beyond the LHB. Study of 
the Earth’s magnetosheath has been carried out through 
in situ measurements of the particle populations and 
magnetic fields. However, there has been no means 
of imaging the entire magnetosheath at the timescales 
necessary for capturing the driving physical processed 
Imaging of the magnetosphere in the soft X-rays pro¬ 
duced by SWCX with either a wide-field imager in Earth 
orbit, or a more traditional instrume nt in solar o rbit can 
cover the entire magnetosphere (Collier et al. 2012). The 
primary issue with such an imager is the effective ca¬ 
dence, which can be determined only in the light of the 
expected flux. Given that the magnetosheath SWCX is 
due primarily to charge-exchange with H, and given that 
that emission is primarily in the \ keV band, then the 


7 Imaging of the magnetosheath with energetic neutral atoms 
ENA) is possible, but only with expo sures of 0.5 to 2 days 
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production factor derived in this paper is essential for 
such calculation. 


1.2. Terminology 

Since there are multiple sites of SWCX, besides comets 
and the atmospheres and surfaces of other planets, it is 
useful to define three different types of SWCX emission 
based on location andphenomenology. 

The Magneto sheatl^j SWCX is due to the solar wind 
plasma interacting with the neutrals, primarily H, in the 
Earth’s extended atmosphere, or exosphere. The terres¬ 
trial magnetic field is an obstacle to the charged parti¬ 
cles that form the solar wind. The magnetopause is the 
surface where the ram pressure of the solar wind is bal¬ 
anced by the pressure of the terrestrial magnetic field. 
Since the solar wind is supersonic and super-Alfvenic, a 
roughly parabolic bow shock forms upstream of the mag¬ 
netopause. The magneto sheath, the region between the 
magnetopause and the bow shock, contains shocked so¬ 
lar wind plasma. Behind the nose of the bow shock, that 
plasma has a density roughly four times higher than that 
in the free-flowing solar wind; as the plasma flows away 
from the nose into the flanks of the magnetosheath, its 
density decreases and its velocity increases. The typical 
stand-off distance of the magnetopause in the direction 
of the Sun is ^ 10 Re from the Earth, but varies with 
the solar wind ram pressure as (nr 2 ) - ^). The magne¬ 
tosheath has a thickness of a few Re that varies with 
solar wind parameters. As the solar wind flux increases, 
the magnetopause is pushed deeper into the exosphere, 
so that the increased solar wind flux interacts with a 
higher neutral density, and the X-ray emissivity increases 
drastically. The continuous compression and expansion 
of the magnetosphere causes the magnetosheath SWCX 
to be strongly time-variable. Because the magnetosheath 
plasma fluxes decrease with distance from the nose of the 
bowshock, the observed magnetosheath SWCX emission 
is also strongly dependent on the observation geometry. 
Simulations of the soft X-ray emissivity in the GSE-XZ[J 
plane are shown in Figure [l] 

The Local Heliospheric SWCX is due to the interaction 
of the solar wind with the neutral interstellar medium 
(ISM) flowing through the solar system. Due to ion¬ 
ization and solar wind pressure, the H component of the 
ISM is mostly excluded from the inner solar system, leav¬ 
ing predominately the He component. The emissivity 
declines with distance from the Sun, due to the R ~ 2 de¬ 
pendence of the solar wind density, so the heliospheric 
emission will be dominated by more local emission. Al¬ 
though the solar wind is highly variable, it is, in the ab¬ 
sence of a strong coronal mass ejection (CME), strongly 
ordered in the Parker spiral pattern. Thus, the local 
heliospheric emission is temporally variable and depen¬ 
dent on the look direction, but not so strongly as the 

8 The term magnetospheric SWCX has also been used (even by 
this author) but is inaccurate. Further, even among heliospheri- 
cists, there is some ambiguity in the meaning and application of 
“magnetosphere”, so readers should be cautious. 

9 All coordinates for the magnetosheath and related regions will 
be given in geocentric solar ecliptic (GSE) coordinates. In this 
right-handed coordinate system, the X-axis is the Earth-Sun line, 
the Z-axis is the ecliptic pole, and the Y-axis is defined by the 
right-handedness. The positive Y-axis is opp osite to t he direc tion 
of the Earth’s motion around the Sun. See Hapgood] (1992) for 
further explication. 
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Figure 1. The magnetosphere at two different solar wind fluxes. Both images are cuts in the GSE-XZ plane. The color scale shows the 
relative soft X-ray emissivity as a function of position. Due to the large dynamic range, the color stretch is problematic, and so has been 
allowed to wrap. The solid white contour traces the region where the field lines are closed, while the large white circle marks a radius 
of 10 R e from the center of the Earth while the small white circle marks a radius of 2.5 Re- The solid green contour traces the region 
removed in the course of removing the plasma not related to the solar wind. In both images the interplanetary magnetic field has B z = —5 
nT (southward pointing) and the time is near the summer solstice. In these images one can see the location of the magnetopause and the 
bowshock, as well as emission due to the interaction of the free-flowing solar wind with the exosphere. Left: n swp v swp = 2.4 x 10 s , a 
typical solar wind flux. Right: n swp v swp = 7.8 X 10 s , a strong solar wind flux. The color bar for the right panel extends an order of 
magnitude higher than the color bar for the left panel. 


magnetospheric emission. A simulation of the soft X- 
ray emissivity in the Earth’s orbital plane is shown in 
Figure [2] 

The Distant Heliospheric SWCX is due to the solar 
wind interacting with both interstellar H and He. The 
pathlengths are long, so the time variability is strongly 
reduced, but should reflect the solar cycle. Due to the 
roughly parabolic shape of the heliosphere, the strength 
must depend upon the look direction. The “ upwind ” 


direction is roughly £ = 5? 1,6 = 19? 6 (Lallement et al. 
2005) which suggests that the ISM flow is due both to the 


solar peculiar motion as well as either radial gas motion 
in the Galaxy or e xpansion from the L oop I superbubble 
(see the review by Frisch et al. 2011). In the following, 
this component will be taken as temporally constant. 

For completeness, it should be noted that there will 
be SWCX due to the heliop ause, bu t this will b e a mi- 
nor component. Using the Koutroumpa et al. (2006) 
model for the distribution of interstellar neutrals m the 
heliosphere and the 0 +7 cross-section therein, one can 
show that the fraction of 0 +7 ions reaching a nominal 
heliopause at ^ 100 au is ^ 50%. If one makes the (un¬ 
likely) assumption that all the remaining ions charge- 
exchange in the heliopause, the ratio of heliopause emis¬ 
sion to the total heliospheric emission is ^ 0.03. Assum¬ 
ing cross-sections of a similar order of magnitude for all 
of the species emitting in the | keV band, the heliopause 
is not a significant contributor to the heliospheric SWCX. 

Recent efforts have attempted to determine which 
of the temporally variable components dominates the 
SWCX “contamination” of XMM-Newton or Chandra 


observations (Kuntz & Snowden 2008; Car ter e t al.|2 
Henley & Shelton 2012; Wargelin et al. 20141). Th< 


2011 


Henley & Shelton ZUIZ; Wargelm et al. ZUI4I). The re¬ 
sults have been equivocal and suggest that both com¬ 


ponents of the emission can be important in different 
situations. The simulations required to determine the 
SWCX production factor bring new insights to this issue 
and will be discussed below. 

1.3. Methodological Overview 

The SWCX emission in a single emission line is given 
by the integral along the line of sight: 


F = 


dCl 

J ( n n n swp v re i(a)fb)—dl (1) 


where n n is the density of neutral targets, n swp is the 
solar wind proton density v re i is the relative velocity of 
the neutrals and the ions, / is the ratio of the density of 
the ion producing the line to the proton density, (a) is 
the velocity-weighted interaction cross-section, and b is 
the fraction of interactions that produce a photon in the 
line of interest. It is necessary to use an averaged inter¬ 
action cross-section since the interaction cross-section is 
velocity dependent. However, the velocity dependence is 
generally small over the range of velocities found in the 
solar wind, so the uncertainty introduced by the need 
to weight by the velocity distribution is generally on the 
order of the uncertainty in the cross-section i tself (see, 
formcample, the catalogue of cross-sections in Bodewits 
2007). The relative velocity of the ions and the neutrals 


is given by 


V rel = (vLp + % 


therm 


( 2 ) 


where v therm 3 kT/nip for the solar wind. 

Semi-empirical models of the distribution of t he neu- 


tral material within the line of sight exist (Fahr 1971 i 
Hod ges|[l994 ), as do semi-empirical models for the solar 
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Figure 2. Critical quantities for the heliospheric model. Top Left: The H distribution in the Earth’s orbital plane in the inner solar 
system. The circle containing data is centered on the Sun and has a radius of 10 au. The downwind direction is marked as are the ecliptic 
longitudes of 0° and 90°. The location of the Earth this simulation step is marked. The color scale is linear and runs from a minimum 
at black/purple (0 cm -3 ) to a maximum at red/white (6.07 X 10 -2 cm -3 ). Top Right: The He distribution, showing the He focussing 
cone. The projection and key is the same as in the previous panel. The color scale is linear and runs from 0 to 5.2 x 10 -2 cm -3 . Bottom 
Left: The proton density in the solar wind. Here the longitude of the ascending node of the solar equator is marked. The color scale 

is logarithmic and runs from 10 —2 cm -3 to 10 3 cm -3 . Bottom Right: The relative soft X-ray emissivity, e = (nn J r nHe) n swpV re l 
assuming that the H and He charge-exchange cross-sections have a ratio of 2:1 The color scale is logarithmic. 


wind (Odstrcil 2003) and measurements of the solar wind 
at the Earth. Thus, one should be able to calculate the 
SWCX emission within a given line of sight at a given 
time, though only to within the rather large uncertainty 
inherent in the models. However, the primary problem 
is that the interaction cross-sections, (cr), and branch¬ 
ing ratios, 6, are poorly measured, particularly in the 
energy range of astrophysical interes t (see t he extensive 


discussion of alternative methods in Smith et al. 


2014). 


However, some important single lines, such as O VII and 
O VIII have been somewhat characterized. The vast 
number of faint blended lines that contribute to the | 
keV band, however, have not. Since, with the notable 


exception o f the Diffuse X-ray Spectrometer ([Sanders 
et al. 2001 [) and X-ray Quantum Calorimeter (McCam- 


mon et al.||2002 ), the | keV band has not been observed 


with sufficiently high spectral resolution to identify line 
complexes, let alone individual lines, we must deal with 
the aggregate of all of these lines. In such a case the | 
keV band emission can be written as 


F = 


where 


noo 

/ ( Tl n Tl SW pV re i^dl 

J 0 


<="2»/A 


( 3 ) 


( 4 ) 
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is the production factor for the energy band, summing 
over all of the different species emitting in the band. It 
is the aim of this paper to determine for the ROSAT | 
keV band. The above expression is often simplified to 

F = sQ (5) 

where 

/*oo 

Q = / {jlnW'swpVrel)dl . ( 6 ) 

J 0 


This expression assumes that the ionization structure is 
constant along the line of sight. Although such an as¬ 
sumption is surely incorrect, we will see the degree to 
which it is immaterial; the magnetosheath is smaller than 
the spatial/temporal scales upon which we can currently 
measure the changes in the ionization structure while the 
heliosphere is so large that the variations disappear when 
integrating along the line of sight. 

Our primary source for | keV data is the ROSAT All- 


Sky Survey (RASS). In the course of constructing the 
maps from the survey, it was realized that there was a 
temporally variable component of the soft X-ray back¬ 
ground whose frequencies ranged from hours to days. 
Longer-term variations could not be characterized with 
the RASS, while shorter term variations were identified 
with othe r phenomena, such as aurorae, and discarded. 


Freyberg (1998) first suggested that these Long-Term 
Enhancements (LTE, |Snowden et al.||1994|) we r e cor- 
related with the solar wind and Cravens et al. (2001) 
demonstrated the correlation more""concretely. Since the 
foreground rate towards the dark Moon was roughly con- 
sistent wit h the LTE rate at the time of the observation 
(Collier et al. 20T4|, it would seem that the LTE were 
due to SWCX in the magnetosheath where the neutral 
species charge-exchanging with the solar wind is predom¬ 
inately H. In this case, one should be able to use an MHD 
model of the magnetosphere for the spatial distribution 
of the solar wind ions in the near Earth environment, a 
model of the H distribution in the exosphere, and mea¬ 
surements of the local solar wind flux to calculate the 
quantity 



'W j swpVrel)dl — Q M - 


(7) 


Here, n en is the density of exospheric neutral H and Qm 
is the Q for the magnetosheath alone. The ratio be¬ 
tween the observed LTE rate and ^Qm should yield 
the production factor, for the ^ keV band for solar 
wind charge-exchange with neutral H. 

There are, of course, a number of complications to be 
considered. The primary complication is the local helio¬ 
spheric SWCX emission, which is due to the interaction 
of the solar wind with (predominately) neutral interstel¬ 
lar He flowing through the solar system. If the LTE flux 
were correlated with the local heliospheric SWCX emis¬ 
sion, we would need to remove the variation due to the 
local heliosphere from the LTE light-curve before deter¬ 
mining the production rate from the magnetosheath. In 
the following, we will use a MHD model of the solar wind 
and a semi-empirical model of the interplanetary distri¬ 
bution of neutral H and He to measure the correlation 
between the local heliospheric SWCX emission and the 
local solar wind flux. We will show that the local helio¬ 


spheric emission provides a negligible contribution to the 
LTE rate. 

To address the issue of the existence of the LHB re¬ 
quires more than the H production factor (d erived here), 
and the He production factor derived by Galeazzi et al. 
(2014), it requires careful understanding of the shape of 
the heliosphere and the extent of emission from the he¬ 
liopause. Both of these issues are beyond the scope of 
this work and will be addressed in the future. 

In the following, §2 describes the LTE data, the so¬ 
lar wind data, the MHD model of the magnetosphere, 
the MHD model of the solar wind for R < 10 au, the 
semi-empirical model of the interplanetary distribution 
of neutral H and He, and the semi-empirical model of 
the Earth’s exosphere. Section 3.1 describes the correla¬ 
tion between the local solar wind flux ( n swp v swp ) and the 
magnetosheath emissivity while §3.2 describes the corre¬ 
lation between the ROSAT LTE and the local solar wind 
flux. Section 3.3 demonstrates that the local heliospheric 
SWCX emission does not contribute significantly to the 
n S wpV S wp -LTE correlation. The discussion in §4 explores 
the implications of the SWCX production factor. That 
section also discusses the implications of the heliospheric 
correlation analysis for observation and mission planning, 
as well as the question of the relative importance of the 
magnetosheath and heliosphere to SWCX “contamina¬ 
tion” . 


2. DATA 
2.1. X-ray Data 

Orbit: ROSAT was launched 1 June 1990 into a nearly 
circular 96-minute orbit with an inclination of 53°. The 
RASS observations began 11 July 1990 and continued, 
with some interruptions, for roughly six months. During 
the All-Sky Survey observations, the spacecraft was spun 
on an axis perpendicular to the optical axis. The spin 
axis was set to be within ~ 13° of the Sun and precessed 
at ~ 4' per orbit, nearly the same rate as the Earth re¬ 
volves about the Sun. The spin period was nearly the 
same as the spacecraft’s orbital period about the Earth. 
Thus, over the course of a single orbit, the spacecraft ob¬ 
served a great circle roughly perpendicular to the Earth- 
Sun line passing through the ecliptic pole^] Over the 
course of the bulk of the survey, the spin axis drifted to 
within ^ 7° of the Sun. 

LTE Isolation: The ROSAT Position Sensitive Pro¬ 
portional Counter (PSPC, the instrument with which the 
RASS was executed) FOV had a diameter of 2°, while 
successive scans were offset by roughly 4'. Thus, each lo¬ 
cation on the sky could be observed for up to 30 consecu¬ 
tive orbits at the ecliptic plane and for many more orbits 
toward the ecliptic poles. Comparison of measurements 
of the same location on the sky from successive orbits al¬ 
lowed the detection and measurement of the time-varying 
background component called the Long-Term Enhance- 

10 Since the orbit precesses by ~ 18' per orbit while the plane of 
the look direction is fixed in GSE coordinates, the relation be¬ 
tween location of the spacecraft and the look direction is con¬ 
tinuously changing. The one constant is that the look direc¬ 
tion is at roughly + (-)GSE-Z when the spacecraft is at its maxi- 
mum(minimum) GSE-Z. Similarly for GSE-Y. GSE-X however is 
more complicated. It should also be noted that the spin vector 
was reversed on a number of occasions to keep the spacecraft from 
scanning the Earth. 
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Figure 3. The 1 keV LTE flux (solid lines) and the solar wind flux ( crosses ) as a function of date. The LTE flux is in units of ROSAT 
count s -1 FOV -1 . The solar wind flux from IMP 8 has been scaled to match the LT E flux. Periods of LTE data for which there is also 
solar wind data are plotted in red. The vertical lines denote the interval considered by |Cravens et ah] ( |2001 ). 

merits (LTE). The details of the identification and mea¬ 


surement of the LTE in the RASS are given in Snowden 
et al.|(|1994| 1995). This time-varying component was re¬ 
moved from the RASS, although some residual LTE con¬ 
tamination remains and can be identified as stripes along 
lines of constant ecliptic longitude. Any non-temporally 
variable component due to the SWCX remains in the 
RASS, and it is the strength of this component that has 
bee n problematic and controversial (compare, for exam¬ 

ple, Lallemen t|2004URobertson et al.|2009| |Koutroumpa| 


et al.||2009[ |WeTsE. et al.||1999p. 

Binning: The | keV band, C band, or R12 band LTE 
light-curve is shown in Figure [3J The light-curve is dis¬ 
continuous because the observations used for the RASS 
were not entirely continuous, mostly due to passage of 
the satellite through the radiation belts and the South 
Atlantic Anomaly. A portion of the | keV L TE light- 
curve, rough ly days 270 to 300, was used by |Craven s 


et al. 


(200l|) to demonstrate that the LTE were corre¬ 
lated with the solar wind flux (as measured by IMP 8 ), 
and thus were due to SWCX. For this work we will use 
the entire light-curve, though the contributions past day 
400 do not contribute significantly. The data were origi¬ 
nally binned into ^ 5 minute intervals with a mean count 
rate of 0.7 count s -1 . No uncertainty was assigned in the 
original measure of the LTE rates and it is clear, given 
the way the values were derived, that the systematic un¬ 
certainties must be larger than the Poisson statistics. We 
have rebinned the data into a uniform 95-minute spac¬ 
ing and have taken the uncertainty of each bin to be 
cr/y/N where a is the RMS variation of the data points 
contributing to an individual bin. Bins containing fewer 
than two data points were discarded. The bin spacing 
was chosen so that each bin contains roughly a single or¬ 
bit of data and that adjacent data points cover nearly the 
same part of the sky. Thus, any systematic effects due to 
the cosmic background will vary by a small amount over 
a day and a half, rather than producing large variations 
between adjacent data points. 

Cusp Removal: For the bulk of each orbit, the line 
of sight passed through the flanks of the magnetosheath 
where the SWCX emission is relatively low. However, 
in some seasons, the line of sight scanned near the mag- 
netospheric cusps, where SWCX should be strongly el¬ 
evated, as shown in Figure [l] The cusps are relatively 


narrow throats of magnetic field lines just poleward of 
the last closed field line on the Earth’s day side where 
the Earth’s magnetic field lines are open and thus they 
are the regions where the solar wind plasma can have 
direct access to the upper atmosphere. The cusps are 
threaded by the field lines undergoing or having recently 
undergone reconnection, and the particle population is 
determined by kinematic effects not included in MHD 
theory. Due to this limitation of MHD models, the den¬ 
sity structures within the cusps are not properly char¬ 
acterized, and so need to be excluded from the analysis. 
Typic ally, the cusp extends 3-4 degrees in mag netic lat¬ 
itude (Zhou et al. 2 0001 Palmroth et a l. 2001) and 2-3 
hours m longitude (|Merka et al. " 2002|) at a given alti¬ 
tude. Due to the curved magnetic field lines threading 
the cusp, a low altitude spacecraft can observe lines of 
sight through the cusp over ^ 15° in latitude, meaning 
that the spacecraft will scan across it in less than five 
minutes. In the RASS processing, short-term enhance¬ 
ments (STE) lasting just a few minutes when the line of 
sight passed near the Earth’s pole were attributed, at the 
time, to auroral emission but may have been due to the 
cusps. These STE periods were removed from processing, 
and thus are not included in the LTE data. 

Removal of the STE periods does not guarantee the 
removal of cusp periods. We have attempted to iden¬ 
tify periods when the line of sight passes through the 
cusp. This is not a trivial matter as the location of the 
cusp varies with diurnal motion, the annual motion, the 
direction of the interplanetary magnetic field, and the 
strength of the solar win d flux. We have used the |Tsy-] 


ganenko & Sitnov (2005) model of the Earth’s magnetic 
field to determine tne"magnetic field geometry along the 
line of sight, that is, whether the magnetic field lines 
through which the lines of sight pass are closed, open, or 
interplanetary. 

We first compared the LTE light-curve with a repre¬ 
sentation of the magnetic field geometry along the line of 
sight and saw no obvious correlation. We further checked 
time intervals for which viewing through the cusp was 
most likely: when the spacecraft was on the day side 
and the line of sight was within 2° of the GSE-Z axis, 
or when we might expect to be under (over) the north- 
ern(southern) cusp and to be looking up(down) through 
it. For each time interval meeting these criteria and for 












































Observations on SWCX 


7 



Figure 4. The | keV LTE flux (solid lines) and the solar wind flux ( crosses ) as a function of date. The LTE flux is in units of ROSAT 
count s -1 FOV -1 . The solar wind flux from IMP 8 has been scaled to match the LT E flux. Per iods of LT E data for which there is also 
solar wind data are plotted in red. The vertical lines denote the interval considered by |Cravens et ah| ( |200 1). 


which there were both LTE and solar wind data, we cal¬ 
culated the difference, A 2 , between the LTE rate at that 
time and the LTE rate at ±10 minutes. Since the scan 
rate is 3?75 minute -1 , the cusp should fall in the center 
bin while the bins 10 minutes away should be free of the 
cusp. The histogram of the values of A 2 is peaked at 
-0.05 count FOV -1 s -1 , has a cr = 0.25 count FOV -1 
s -1 . There is a slight positive asymmetry, with an ex¬ 
cess of events with 0.5 to 1.0 FOV -1 s -1 , and with a 
values extending to 1.3 count FOV -1 s -1 . Time inter¬ 
vals with particularly high A 2 did not show any partic¬ 
ular distribution in spacecraft coordinates. Since STEs 
were removed from the LTE data and our search for en¬ 
hancements at times likely to contain cusp emission did 
not reveal systematic increases in emission, we can con¬ 
fidently state that the residual contribution by the cusp 
to the LTE light-curve is small. 

The | keV Band: The M band, or R45 band LTE 
light-curve is shown in Figure [ 4 J The | keV band data 
were treated in the same manner as the C-band data. 
Comparison to Figure [3] shows that the LTE in the two 
bands are rather different. The differences are due to 
the SWCX spectrum. Although the details of the SWCX 
spectrum are poorly understood, the overall shape of the 
spectrum is set by the ion abundances in the solar wind 
and the location of their principal recombination lines. 
The | keV band is dominated by the O VII and O VIII 
lines with occasional significant contributions from Ne IX 
and Mg XI (see, for example, Carter & Sembay 2008). 
Rather than a few strong lines, the ^ keV band contains 
a large number of faint lines from many different species. 
Thus, the | keV band light-curve strongly depends upon 


the relative abundance of oxygen (compared to protons) 
and the relative ionization states of that oxygen. Con¬ 
versely, since the | keV band contains such a multitude 
of lines, none of which dominates the band, the abun¬ 
dance and ionization variations average out. 


2.2. Solar Wind Data 

Solar wind data for the period of the RASS are some¬ 
what sparse for two reasons. First, the only satellite 
monitoring the solar wind at that time was IMP 8. Since 
IMP 8 was launched to study not only the solar wind, 
but the magnetospheric boundary and the magnetotail, 
it was exposed to the free-flowing solar wind for only 


a portion of its orbit. IMP 8 did not have a means of 
recording data, and telemetry coverage was only 65% to 
80% in this period. Thus, the solar wind parameters 
are measured for only a portion of the period for which 
we have LTE data. The solar wind data themselves were 
drawn from the NASA OMNI compilatiorj 11 ] as that data 
set has been cleaned of inappropriate data periods and 
other data artifacts, as well as time shifted in a stan¬ 
dardized manner. From these data we derive the “local 
solar wind flux”, that is, the solar wind flux, the prod¬ 
uct of the proton density and the proton speed, in the 
neighborhood of the Earth. 


2.3. Magnetospheric Simulations 

We represent the ion density of the magnetosheath 
with simulations using the BATS-R-US (Block-Adaptive- 
Tree Solar Wind Roe-Type Upw ind Scheme TM code 
( Toth et al.|2005 Powell et al. |1999 ). BATS-R-UB solves 
the three-dimensional 1V1HD equations for the region sur¬ 
rounding the Earth, simulating the response of the mag¬ 
netosphere to the temporally variable solar wind. BATS- 
R-US is one of several standard MHD models of the mag¬ 
netosphere in common use. Validation of these models 
by comparison with in situ spacecraft measure ments is 
an ongoing effort among many institutions (e.g. Raeder 
2003, among others) and a new intercomparison of mod¬ 
els is in progress (Collado-Vega & Sibeck, in preparation, 
2015). 

It should be noted that the magnetosheath is sensitive 
not only to the solar wind flux and the IMF but also 
to the recent history of those parameters. This depen¬ 
dence is due, in part, to the time required for a solar 
wind impulse to move from the nose of the bowshock 
past the Earth (^ 3 minutes) and down the tail. The 
dependence is also due to the finite time for the magne¬ 
tosphere to respond to changes in the solar wind. Thus 
MHD mo dels are s uperior to analytic models, such as 
that of Spreiter et al. (19 66|), which have been used for 
previous SWC X studies ([Carter et al. _ 20111 Knntz fe 
Snowden]|2008 ). In comparison to Spreiter et al. ( |l966| ) 


11 ftp://spdf.gsfc.nasa.gov/pub/data/omni/high_res_omni/ 

12 Simulation results have been provided by the Community Co¬ 
ordinated Modeling Center at Goddard Space Flight Center. The 
BATS—R—US model was developed by Dr. Tamas Gombosi at the 
Center for Space Environment Modeling, University of Michigan. 
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Table 1 

BATSRUS Runs 


Name 

Version 

Dipole 

Start 

Step 

Number 

Label 

of Run a 


Tilt 

Date 

Size 

of Steps 




(deg) 


(min) 




Northern hemisphere winter solstice 

SWPC_SWMF_052811_2 

V20110131 

-32.50 

2006/12/14 

1 

2190 

O 

Marc_Kornbleuth_120513_3 

V20110131 

-31.60 

2002/12/19 

10 

37 

X 

Marc_Kornbleuth_l 11213_3 

V20110131 

-26.26 

2004/11/06 

10 

6 

□ 

Ankush_Bhaskar_050612_l 

v20110131 

-25.81 

2012/01/08 

5 

139 

+ 

Equinox 

Yaireska_Collado-Vega_l 12812.1 

V20110131 

-2.50 

1991/03/24 

1 

601 

O 

SWADESH.PATRA.040412.2 

v8.01 

-2.06 

2002/03/24 

4 

510 

n 

Lur_Zizare_032112_l 

v8.01 

-0.45 

2012/03/12 

4 

511 

X 

tao_huang_082814_6 

V20140611 

-0.17 

2008/02/23 

4 

76 

A 

Ilja_Honkonen_073014_2 

v20130129 

2.15 

2005/08/22 

4 

1258 

+ 

Northern hemisphere summer solstice 

Yuni_Lee_012610_3 

v8.01 

19.00 

2001/06/01 

5 

132 

A 

SWPC_SWMF_060411_6 

v20110131 

29.50 

2005/05/14 

12 

388 


Steven_Snowden_011212.1 

v8.01 

28.66 

2006/06/04 

5 

337 

n 

Yaireska_Colladovega_091112.1 

V20110131 

28.88 

2007/05/20 

1 

50 

o 

Chigomezyo_Ngwira_022014_l 

V20130129 

29.83 

2013/06/28 

1 

1200 

X 

Brian_Walsh_030413_l 

v20130129 

30.04 

2000/07/15 

1 

338 

+ 


a The name of the run as it appears on the CCMC run-on-request website. 


we find that BATS-R-US produces consistently different 
positions for the magnetopause and the bowshock (see 
§4.2). BATS-R-US also models the asymmetry of the 
magnetosheath. 

Each BATS-R-US simulation is sampled at a predeter¬ 
mined cadence, usually every few minutes. Each sample 
produces proton density, speed, and temperature values 
for each point in a three-dimensional grid whose spacing 
is optimized to provide high resolution of strong gradi¬ 
ents. One disadvantage of the BATS-R-US code is that 
it does not distinguish between protons originating in the 
solar wind, and those much lower energy protons origi¬ 
nating in the plasmasphere closer to the Earth. However, 
the amount of solar wind pl asma wit hin the magneto¬ 
sphere is relatively low (Christon et al. 1994). Thus, one 
can isolate the solar wind protons by removing regions 
that lie within the Earth’s closed field lines. This method 
has difficulties around the cusps, but as the cusps are 
problematic with MHD models, and have been removed 
from the LTE light-curves, this issue is not significant for 
this analysis. (However, see §3.1.) 

Running BATS-R-US for the roughly 200 days of the 
RASS would be a significant effort, particularly with¬ 
out continuous solar wind data. At the LTE sample 
rate of five minutes, we would accumulate 57600 sam¬ 
ples, at roughly 75Mb each. Even sampling once per 
orbit would produce 3200 samples. While such sampling 
would be ideal, it is beyond the scope of this work. We 
have, however, accumulated a series of BATS-R-US sim¬ 
ulations from other projects that cover a variety of solar 
wind conditions. The runs and their relevant parameters 
are listed in Table [l] The validity of using representa¬ 
tive simulations rather than dedicated simulations of the 
RASS epochs will become apparent in the analysis. 


2.4. The Model of the Exospheric Neutral Distribution 


The neutral densi ty inthe magnetosphere is taken from 
the model of Hodges (1994). This static model was orig¬ 
inally calculated for a variety of values of insolation as 
measured by flO.Tpd at both solstice and equinox. For 
any particular set of conditions and dates, we interpo¬ 
late among the available model states. The model is valid 
from 1.05 R#to 9.75 R#; at larger distances the model is 
extrapolated as R~ 3 . Comparison of the Hodges model 
calculations wit h Lyman a column b rightnesses for the 
antisolar point (0stgaard et al. 2003) suggests that the 
Hodges model is compatible with the limited observa¬ 
tions. 


2.5. Heliospheric Simulations 

We represent the ion density in the heliosphere with 
output from the ENLIlj^j code (Odstrcil 2003). ENLIL 
solves the three-dimensional M H DWqu at ions lor the in¬ 
ner heliosphere. The inner boundary conditions of the 
solar wind density, speed, and temperature are derived 
from solar magnetograms. The inner boundary is taken 
to be 21.5 R 0 where the solar wind becomes supersonic. 

The ENLIL simulation used for this work was specially 
requested of the CCMC. The first sampled time step of 
the simulation is 1990-09-28 12:18 UT. The simulation 
was sampled with a time step of 90 minutes (roughly 
the ROSAT orbital period) for a total of 529 samples, or 
slightly over a mo nth of simulation, cor responding to the 
period studied by Cravens et al. (2001). The simulations 
were run over a standard region; 21.5 R 0 < R < 10 


13 The 10.7 cm (2800 MHz) solar flux as measured by the Pen¬ 
ticton observatory. 

14 Simulation results have been provided by the Community 
Coordinated Modeling Center at Goddard Space Flight Center. 
ENLIL was developed by Dusan Odstrcil at the University of Col¬ 
orado at Boulder. ENLIL is not an acronym, despite its nonstan¬ 
dard capitalization, but the name of the Sumerian god of winds 
and storms. 
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au and 60° < 0 < 120° where 0 is the angle from the 
north solar rotational pole. The restriction to within 30° 
of the solar equator will not be of grea t concern for the 


correlation analysis to be done in {3.3 


2.6. The Model of the Heliospheric Neutral Distribution 

The neutral density in the heliosp here is rep resented 
by the model developed by Koutroumpa et al. (2006) 


based on the classical hot model calculations of 


ment et al. (11985), Lallement et al. (1985), and 


Lallc- 


The model was calculated for 


Lalle- 

solar 


ment et al.[ 

maximum conditions such as those that would have been 
experienced during the ROSAT All-Sky Survey obser¬ 
vations. The ratio of photon pressure to gravitational 
force on H, /i, is set to 1.46, while the total equato¬ 
rial ionization rate (including CX with protons, pho- 
toionization and maybe electron impact ) derived from 
the SWAN data (Quemerais et al. 2006) is 6.54 x 10 -7 
0-1 The H ISM para meters used for construct ing these 


models are taken from |Lallement et al. (2005): upwind 
(A,/3) = (252?3,8?5), T = 13,000k,A;o = 21.0 km 
s -1 , and no = 0.1 cm -3 . The equatorial electr on im¬ 
pact factor for He is given by the rates in |Rucinski fe 
Fahr| fl 19891) mul t iplied by 2.05, the factor estimated by 
LallementTet al. (2004) to account for the increase dur¬ 
ing s olar maximum. The H e IS M parameters are taken 
from Mobius et al. ([2004) and Witte (2004): upwind 
(A,/3) = (254?7, 5?2), T = 6300K7^o = 26.3 km s -1 , 
and no = 0.015 cm -3 . The H and He contributions are 
calculated separately as the ion-neutral interaction cross- 
sections are different for the two species. 


3. ANALYSIS 

We will first determine the relation between the local 
solar wind flux, n swp v swp and the Qm = f ^en^swp^reidl 
for the magnetosheath. It is expected that this relation 
should be non-linear; as the solar wind flux increases, 
the bowshock is pushed closer to the Earth, where the 
neutral density is higher. We must determine what part 
of this relation is comparable to the relation between 
n SW pV S wp and the | keV LTE rate. By combining the 
two relations, we then get a direct relation between Qm 
and the \ keV LTE rate, and thus the ROSAT \ keV 
band flux. 


3.1. The Qm - Solar Wind Flux Relation 

In order to understand the magnetosheath SWCX flux 
that ROSAT ought to have observed during the All-Sky 
Survey, we sampled a number of magnetospheric simu¬ 
lations with lines of sight typical of the RASS. It is im¬ 
portant to bear in mind that the LTE data with which 
we are to compare the results of simulations are binned 
to intervals of roughly one orbit. Thus we wish to ex¬ 
tract from the simulations the equivalent to what ROSAT 
would have observed over a single great circle scan. We 
begin with a “reference” scan, a great circle, centered 
on the Earth, passing through the ecliptic poles, and ro¬ 
tated 12° about the GSE-Z axis from the GSE-YZ plane. 
For each time step of each simulation we calculated in¬ 
tegrated emissivity, Qm, for 360 lines of sight uniformly 
distributed over this great circle. From those we then 
calculated the mean integrated emissivity Qm for the 


scan. Since this reference scan passes through the eclip¬ 
tic poles (the GSE-Z axis) it avoids the cusps except for 
times aro und the solstices. Once we understand the be¬ 
havior of Qm for the reference scan, we will explore the 
effects on Qm of using more realistic scans. 

As noted in §2.3, the BATS-R-US simulations do not 
distinguish between solar wind protons and the cold plas- 
maspheric protons, while the SWCX emission can be as¬ 
sociated with only solar wind protons. We thus need to 
remove the plasmasphere from the simulation before cal¬ 
culating Qm • Since these scans avoid the cusps, the pro¬ 
ton density for any direction on a scan is peaked at the 
Earth, declines with increasing radius to some relative 
minimum, and then shows another peak at the magne¬ 
tosheath. The plasmasphere can be crudely but rather 
efficiently defined as the region inside the relative mini¬ 
mum in proton density, though the plasmasphere by no 
means fills that region. Figure [5] shows the proton den¬ 
sity for a GSE YZ slice through one step of a simulation, 
as well as the boundary of the “plasmaspheric” region to 
be removed. This cleaning algorithm works well when 
the cusp does not enter the scan. 

In the top two panels of Figure [6j each plot point repre¬ 
sents the Qm for an individual time step of a simulation 
listed in Table |TJ The points have been color-coded by 
season and, where possible, symbol-coded by simulation. 
The top panel shows the bulk of the parameter space 
covered by the simulations, while the middle panel cov¬ 
ers only the values of n swp v swp found in the LTE data 
periods. The bottom panel shows the values binned into 
An SW p v swp = 0-2 x 10 8 cm -2 s -1 wide bins. The colored 
symbols are points from individual simulations while the 
black points are summed over all the simulations. The 
standard deviations are shown for bins containing more 
that five values. 

As one might expect, as the solar wind flux increases, 
the value of Qm increases non-linearly. The relation be¬ 
tween Q and n swp v swp becomes steeper as n swp v swp in¬ 
creases. Below n swp v swp = 6 x 10 8 cm -2 s -1 the re¬ 
lation is relatively linear, or at least any non-linearity 
is obscured by the intrinsic variation. The bottom 
panel of Figure [6] shows a linear fit to the values with 
n swp v swp < 6 x 10 8 cm -2 s -1 and that fit seems ade¬ 
quate. The values do not deviate significantly from the 
linear fit until n swp v swp 10 x 10 8 cm * 2 s 1 . There 
is a significant dispersion in Qm at any given value of 
nswpVswp for a single simulation, and clear systematic 
differences between simulations. At n swp v swp > 10 x 10 8 
cm -2 s -1 there are fewer simulation points and the be¬ 
havior of the relation becomes quite uncertain. 

The first source of variation in Qm at a given n swp v swp 
is the inherent difficulty of defining the solar wind pres¬ 
sure acting upon the magnetosheath. The solar wind 
takes roughly three minutes to move from the nose 
of the bowshock to a plane passing through the cen¬ 
ter of the Earth. Each portion of the magnetosphere 
responds to the ambient pressure at its location, and 
the response propagates rapidly throughout the magne¬ 
tosheath. Thus, setting a single characteristic n swp v swp 
for a time step is inherently inaccurate. There is, as 
well, the difficulty of propagating strongly variable solar 
wind conditions from their point of measurement to the 
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Figure 5. Left: A cut through the magnetosphere in the GSE YZ plane. The colorscale traces the proton density and shows a peak 
around the Earth (the plasmasphere) and a peak at 20-25 Re (the magnetosheath). The white circle is centered on the Earth and has a 
radius of 18 R^;. The black line shows the boundary of the region excluded by the plasmasphere removal algorithm. The north ecliptic pole 
is up. Right: A cut through the magnetosphere in the reference scan plane (Y'Z) which is rotated 12°about the GSE-Z axis from the 
GSE-YZ plane. The colorscale traces the SWCX emissivity. The maximum emissivity is 2.3 X 10 9 cm -5 s -1 . The black circle is centered 
on the Earth and has a radius of 18 R#. The contours trace the logarithm of the neutral density. The black region at the center is the 
plasmasphere as defined by the plasmasphere removal algorithm. 


(moving) magnetopause or bowshock. Thus, there is an 
intrinsic error in setting n swp v swp , particularly when the 
solar wind conditions are changing rapidly. Since high 
n SW p v swp periods tend to be periods with strong, rapid 
variation in n swp v swp , the variation in Qm increases with 
nswpVswp • The simulations with slower changes in the so¬ 
lar wind show significantly less scatter. 

A second source of variation is the asymmetry of 
the magnetosheath, which depends not so much on the 
n SW p v swp of the solar wind, but upon the variation in 
the direction of the solar wind momentum vector and 
the direction of the interplanetary magnetic field. The 
asymmetry can be seen in Figure |5j Nor is the neutral 
density spherically symmetric; rather it is squashed at 
the poles, which leads to a rather asymmetric distribu¬ 
tion of emission. Since the scan path is rotated by 12° 
around the GSE-Z axis with respect to the YZ plane, it 
samples the magnetosheath at much closer distances on 
one side than it does on the other, which leads to even 
more extreme asymmetries. The asymmetry in the emis¬ 
sion for a typical scan, after removal of the plasmasphere, 
is shown in the right-hand panel of Figure [5] The aver¬ 
age of the integrated magnetospheric emissivity over a 
great circle scan, Qm , depends on the direction and am¬ 
plitude of the rotation of the scan plane with respect to 
the underlying asymmetry. Thus, for any given value of 
nswpVswp, one would expect a variation in Qm due to the 
variation in the asymmetry. Using a representative set 
of simulations, reversing the sense of the rotation of the 
scan path with respect to the GSE-YZ plane changed the 
Qm for individual time steps from -20% to +10%, with 
a mea n change of -7%. This variation might be reduced 
if the (Hod ges| [l994) model overestimates the equatorial 
bulge of the exosphere. 

A third source of variation is seasonal variation. As the 


Earth nods, the cusp latitude changes, bringing the cusps 
closer or further from the scan path. However, as indi¬ 
cated in Figure [6j this variation seems to be small com¬ 
pared to the intrinsic differences between simulations. 

The combined effect of all of these variations produce 
that scatter seen in Figure [6j We have calculated the 
fractional variance as a function of n swp v swp . Using val¬ 
ues for all of the simulations, we calculated the mean Q 
for A n swp v swp = 0.2 x 10 8 cm -2 s -1 wide bins (the black 
boxes in the bottom panel of Figure [6]). The standard 
deviation of the values in each bin is snown as the black 
histogram. For n swp v swp < 8 x 10 8 cm -2 , the standard 
deviation in Q, erg, is roughly 40% of the mean Q and is 
relatively independent of n swp v swp . 

By fitting the data in the last panel of Figure [6] we find 
that 


10 8 


= (0.037 ± 0.189) + (20.68 ± 0.41) 


Qm 

10 20 


__ ( 8 ) 
where n swp v swp is in cm 2 s 1 and Qm is in cm 4 s 1 . 
The correlation is fitted from the values binned over all 
simulations (the black boxes in the bottom panel of Fig¬ 
ure |6| since the individual values derived from each time- 
step of each simulation have no intrinsic uncertainties. 
The measured standard deviation of each point was used 
as the uncertainties in the fit. We fit only the values be¬ 
low n swp v swp < 6 x 10 8 cm -2 s -1 where the relation is 


most clearly linear. We have used a fitting method (Akri- 


tas & Bershady 1996) that minimizes the uncertainty 
perpendicular to the fit, and determines the uncertain¬ 
ties in the fit parameters using a bootstrap technique. 
Note that the uncertainty in the intercept is significantly 
larger than the intercept itself, thus the intercept is not 
significantly different from zero. 
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Figure 7. The ratio of the Qm for the reference scan to the Qm 
from four test scans, as a function solar wind flux. The test scans 
are color coded by the orbit location at which the scan is in the 
+GSE-Z direction: blue: noon, green: dusk, yellow: midnight, 
and red: dawn. The behavior is expected given that the day-side 
emission is stronger than night-side emission. 

Any intrinsic scatter in the n swp v swp -QM relation 
translates directly into a scatter in the n swp v swp - LTE re¬ 
lation. Besides the scatter explored above, there is also 
scatter due to the difference between the geometry of the 
reference scan and real ROSAT scans. Over the course of 
an orbit, the spacecraft will move from ~ 1.09 Re on the 
day side of the Earth to ~ 1.09 R# on the night side of the 
Earth; the scan path moves correspondingly with respect 
to the reference scan. As the emission increases towards 
the sub-solar point of the magnetosheath, the more sun¬ 
ward the scan, the higher the emission. A real ROSAT 
scan will sample the magnetosheath on both sides of the 
reference scan, but since the emission is not linear with 
distance from the sub-solar point, the effects do not can¬ 
cel. 

To measure the difference between realistic scans and 
the reference scan, we constructed four test scans for a 
single (rather unrealistic) test orbit lying in the GSE- 
XY plane. Each test scan was constructed such that 
the line of sight is towards the north ecliptic pole when 
the spacecraft is at one of the cardinal directions, dawn, 
dusk, noon, or midnight. These four test scans sample, 
however sparsely, the possible range of orbit attitudes. 
We have compared the Qm for these test scans to the 
Qm for the reference scan. The result for a simulation 
near summer solstice when the solar wind flux spanned 
the bulk of the expected range is shown in Figure [7| the 
bulk of the differences between the real scans and the 
reference scan lie between -10% and +20%. 


Figure 6. The n swp v S wp-QM relation for the simulations in Ta- 
ble[l] They are color-coded by season: blue: Northern hemisphere 
winter solstice, green: equinox, and red: summer solstice. Table 
^provides a key to symbols where used. Top: The bulk of the 
parameter space covered by the simulations. Each point is a single 
time step of one simulation. Middle: A detail of the previous 
panel showing the parameter space covered by the LTE data. Bot¬ 
tom: The data binned by n swp v swp . The colored data points are 
the binned data from individual simulations. The standard devia¬ 
tion of the data points in each bin is shown if there were more than 
five data points in the bin. The black points are the binned data 
from all the simulations. The black histogram shows the standard 
deviation for the black points. 


3.2. The LTE - Solar Wind Flux Relation 
The correlation of the | keV band LTE flux with the 
local solar wind flux is shown in Figure [8] (Left) where, 
as in the light-curve, the data have been binned in 95 
minute bins. The uncertainties shown are the RMS vari¬ 
ations of the data being binned. As can be seen from 
the histogram of the solar wind data for the periods for 
which there is LTE data (Figure [8]) the bulk of the data 
is for n swp v swp < 5 x 10 8 cm -2 s -1 . Above that value 
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Figure 8. Left: The 1 keV LTE flux as a function of the local solar wind flux. The solid line is the best fit of (LTE| n swp v swp ), the 
dashed line is the best fit minimizing the distance orthogonal to the fitted line, and the dotted lines show the relations with 0.5 and 1.5 
times the best fit slope. Right: Histogram of the solar wind flux (black) and solar wind speed (red) values for the time periods with LTE 
data. The dotted lines are binned by orbit while the solid lines are for individual five minute intervals. 


of the local solar wind flux the data become increasingly 
sparse, tailing off by 10 9 cm -2 s -1 . The bulk of the LTE 
observations were made when n swp v swp was well within 


the range for which the ti swp v swp -Qm relation is linear, 
and only the most extreme values fall in the regime where 
the linearity is dubious. We have fitted the correlation 
for < 6 x 10 8 cm" 2 s -1 us ing the bootstrap 


method of |Akritas fc Bershady| (1996). We find 


LTE = (0.0012 ± 0.0353) + (0.291 ± 0.014) 


TtswpV swp 

lb® . 
( 9 ) 


where LTE is in count s 1 FOV 1 and n swp v swp is in 

— 2—1 

cm s . 

The scatter in the n swp v swp -LTE relation is not unex¬ 


pected. The measured scatter in the ti swp v swp -Qm rela¬ 
tion is roughly 40%. The difference between the reference 
scan path and a “real” scan path is roughly 10%. As¬ 
suming that these variations are independent, that sug¬ 
gests a roughly 50% variation. The blue lines in Figure [8] 
show this variation. It encompasses the bulk of the data 
points, suggesting that there is no significant extra source 
of variation in the n swp v swp -LTE relation. 

We further note the lack of a significant intercept in 
the n swp v swp -LTE relation. This result is actually rather 
unexpected. The measurement of the LTE rate was 
done without reference to the solar wind flux as that 
relation was unknown at the time. One might expect 
the LTE rate to be systematically underestimated or 
over-estimated and indeed, there are a number of bright 
streaks in the LTE-subtracted RASS that suggest that 
the LTE was underestimated in this process. The lack 
of an intercept in the n swp v swp -LTE relation suggests 
that the estimation of the LTE, while imperfect for some 
scans, was not systematically low or high. 

Therefore, we find that 


C R 12 oc (6.02±0.31) x 10 -20 (count FOV -1 cm 4 )Q (10) 



Figure 9. The | keV LTE flux as a function of the local solar 
wind flux. The dashed line is the best fit. This fit was used for 
scaling the solar wind flux to the LTE data in Figure^] 


where Cr \2 is the ROSAT | keV band (or R12 band) 
count rate in count s -1 FOV -1 and Q has its usual units 
of cm -4 s -1 . The vignetting of the ROSAT PSPC pro- 
dn ces an ef f ective area for the FOV of 1.56 square degrees 
(Snowden et al.J 2015|). (This value can be found directly 
from the ROSAT detect or map for the R12 band.) There¬ 
fore we find that 

^ = (3.86 =b 0.20) x 10 -2 °count degree -2 cm 4 (11) 

The situation for the | keV band is rather different. 
Figure EH shows the correlation of the | keV LTE flux 
and thMocal solar wind flux. Although there are more 
low LTE flux data points at low local solar wind values, 
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there is little correlation between the two values, and the 
Pearson correlation coefficient is only 0.45. In compar¬ 
ison, the Pearson correlation coefficient between the | 
keV LTE flux and the solar wind flux is 0.74. The Pear¬ 
son correlation coefficient between the | keV and | keV 
band LTEs is only 0.14, so the | keV band LTE emission 
is better correlated with the local solar wind flux than it 
is with the LTE emission in the | keV band. The over¬ 
all LTE flux in the | keV band is a factor of seven to 
eight lower than in the \ keV band, so the uncertainties 
must consequently be larger. However, we note that the 
RMS/mean for the | keV band is similar to that of the 
- keV band (33.06 versus 34.67) so the decrease in the 
correlation coefficient is not due to an increase in an in¬ 
trinsic scatter in the | keV band. The lack of correlation 
in the | keV band, compared to the ^ keV band, sug¬ 
gests that the oxygen lines that dominate the | keV are 
less well correlated with the solar wind flux than the ag¬ 
gregate of the many lines that produce the | keV band. 
This is not unreasonable since the oxygen lines will be 
sensitive to the variation in the oxygen abundance and 
ionization fraction in the solar wind, while the aggregate 
of lines in the | keV lines should be less sensitive. 


3.3. The Qh - Solar Wind Flux Relation 

Combining the n swp v re i~ LTE and n swp v re i~Q relation 
to get the production factor is only valid if the only con¬ 
tribution to the correlation between LTE emission and 
the local solar wind is the magnetosheath emission. So 
long as it is not correlated with the local solar wind flux, 
the heliospheric SWCX will introduce only noise to the 
relation. Thus, we must investigate the possible correla¬ 
tion of the heliospheric SWCX with the local solar wind 
flux. 

The heliospheric emission along the line of sight is 
Qh = f ( n i P H + Tn ip He)n S w P Vreidl, where n ipH is the 
interplanetary H density, n ip H e is the interplanetary He 
density, and T is the ratio of the He charge-exchange 
cross-section to the H charge-exchange cross-section. 
Since He cross-sections are usually only a factor of a 
few smaller t han those for H (see the compilation in 
Koutroumpa et al.| 2006), we have set T = 0.5. The solar 
wind proton density was extracted from the ENLIL mod¬ 
els while the density of the interplanetary neu trals was 
taken from the model of Koutroumpa et al. (2006). Since 
the ENLIL model was calculated only to 10 au from the 
Sun within 30° of the solar equator, we can calculate a 
line of sight from Earth to 9 au for only a limited number 
of directions. The distance of 9 au is sufficient for deter¬ 
mining the correlation between the heliospheric SWCX 
and the local solar wind flux: if we assume a heliopause 
distance of 100 au, then 70% of the total observed he¬ 
liospheric emission occurs within the first 9 au. There 
is some variation between the upwind and downwind di¬ 
rections, but we have calculated this as an average over 
angle from the upwind direction. 

For the region for which we can measure the helio¬ 
spheric properties over lines of sight 9 au long, we have 
calculated Qh at 5° spacing in both ecliptic latitude and 
ecliptic longitude for each ENLIL time step. For each line 
of sight we calculated the Pearson correlation coefficient 
between the Qh and the local n swp v swp , as recorded in 


the model at the location of the Earth. Figure [To] shows 
the correlation coefficient as a function of ecliptic lat¬ 
itude and angle from the Sun along the ecliptic. The 
pattern shown in this map is reasonably independent of 
the details of the recent solar wind history and the time 
of year. 

The strongest correlations are found between 100° and 
130° from the Sun, measuring in a positive sense around 
the ecliptic, or 50°to 80°from the anti-sun as seen in 
Figure [lO] This is the region where one is looking along 
the Parker spiral beyond the Earth’s orbit. Similarly, 
angles further than —130° from the anti-sun have high 
correlations because one is looking down the Parker spiral 
within the Earth’s orbit. The curvature of the spiral is 
stronger within the Earth’s orbit than without, so there 
is a shorter distance at which the line of sight is tangent 
to the spiral. Only a small part of the Parker spiral region 
is accessible to either ROSAT or XMM-Newton , so one 
should not expect a strong correlation between the local 
solar wind and the heliospheric SWCX. This statement 
excludes, of course, strong CMEs. 

The restriction to within 30°of the solar equator is not 
of great concern for the co rrelation analysis . As c an be 
seen from the diagrams in McComas et al. ([2003), dur¬ 
ing solar minimum, the solar wind between the equator 
and ^ 30° is characterized by a dense, slow “equato¬ 
rial” flow, while at higher latitudes, the there is a faster, 
more tenuous “polar” flow. As can be seen in Figure [TT] 
there are coherent structures extending out of the plane. 
However, unlike the Parker spiral, these structures have 
a wide range of tilts, so the direction in which one looks 
along the structure is not fixed. Further, for a line of 
sight towards the poles, the equatorial flow accounts for 
only ~ 30% of the heliospheric SWCX flux, assuming an 
emissivity that declines as r -2 as measured from the Sun. 
Thus, in general, the correlation between the local solar 
wind and the SWCX observed towards the ecliptic poles 
should be weak, though SWCX emission will be corre¬ 
lated with the local solar wind for some limited direc¬ 
tions at some times. During the solar maximum, which 
is more relevant for the RASS observations, the differ¬ 
entiation between equatorial and polar flows is missing, 
and the solar wind speeds and densities are strongly vari¬ 
able in both regions. In this case we again do not expect 
correlation between the locally measured solar wind and 
the solar wind at the higher solar latitudes not covered 
by the ENLIL model. 

4. DISCUSSION 
4.1. On the production factor 

In §3.1 and §3.2 we fitted the n swp v swp -LTE relation 
and the n swp v swp -QM relation, and used the slopes of 
those relations to determine the ROSAT ~ keV count 
rate as a function of Qm, and thus the | keV produc¬ 
tion factor. Further, we have shown that the scatter in 
the n swp v swp -LTE relation is comparable to that in the 
n' m. relation. 


Robertson et al. (2009) calculated ROSAT production 
factors for a “slow” solar wind from the available atomic 


data. The slow solar wind conditions (see von Steiger 
et al. 2010 for definition) are appropriate for the equa- 
tonaf now at solar minimum. Although the solar wind 
experienced by the Earth at solar maximum is a complex 
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Angle from Anti—Sun 


Figure 10. The correlation of the heliospheric SWCX flux and the local solar wind flux as a function of look direction. The longitudinal 
coordinate is the angle from the anti-Sun while the latitude is the ecliptic latitude. The color scale and the contours are the correlation 
coefficient between the local solar wind flux (as measured at the Earth in the ENLIL model) and (uh + XnH e ) n swpV r el integrated along 
the line of sight. The two sets of vertical lines mark the locations that can be observed with ROSAT ±15° from the perpendicular to the 
Sun; XMM-Newton is restricted only to ±20° from the perpendicular to the Sun. The effect of the Parker spiral is apparent at 50°to 70° as 
well as between -130°and -150°. This map is the mean of four maps made from the same ENLIL simulation at 90° intervals around the 
Earth’s orbit. Each individual map is very similar to the mean shown here. 



Figure 11. The relative solar wind proton density from a vertical 
cut through a time step of the ENLIL simulation. The circle has 
a radius of 10 au. The location of the Earth is marked. 


mixture of slow and fast solar wind, the distribution of 
velocities during the time of the RASS is consistent with 
the slow wind, and it is common practice to approxi¬ 
mate the equatorial solar wind, even at solar maximum, 
with the slow solar wind ion/abnndanc e ratios. The cal- 
culated \ keV production factor from 


Robertson et al. 


(2009) for charge-exchange with H is 23.64 x 10 count 
arcmin -2 cm 4 , (0.851 x 10 -20 count degree -2 cm 4 ) which 
is a factor of 4.5 lower than that derived here. 

This disagreement is not altogether surprising given 
the uncertainty in the atomic data. A comparison, for 
example, of the soft X -ray SWCX spectrum shown in 


Robertson et al. (2009) with that of Koutroumpa et al. 
foof shows strong differences despite a reliance on 
similar atomic data. A conversion of the IKontroumpaj 


et al.| (2009) spectrum into a production factor yields 
2.13 x 10 20 count degree 2 cm 4 , which is a factor of 
only 1.79 lower than that empirically obtained from the 
LTE data. Addition of further lines to the theoretical 


spectrum will decrease this difference. 

We do not find any combination of systematic uncer¬ 
tainties that could push the measured production factor, 
<N significantly lower. Had there been a heliospheric con¬ 
tribution to the n swp v swp -LTE relation, then the derived 
would have been an upper limit. However, as shown 
in §3.3, the heliospheric emission does not contribute to 
the n swp v swp -LTE correlation. Conversely, if the plasma- 
spheric contribution in the BATS-R-US simulations were 
not completely removed, then the derived ^ would be a 
lower limit. As can be seen in Figure [5j the region re¬ 
moved for the plasmasphere, while adequate for that pur¬ 
pose, does not extend all the way to the magnetopause. 
An alternate cleaning which removes both the plasmas¬ 
phere and the bulk of the low density region interior to 
the heliopause does not significantly change the calcu¬ 
lated Qm, though some individual lines of sight changed 
significantly. Thus, the plasmasphere removal method 
does not significantly change the derived production fac¬ 
tor. Similarly, if the Hodges (1994) mod el overestimates 
the e xospheric density, as suggested by |0stgaard et ah 


(2003), then the derived ^ would again be a lower limit. 

The derived production factor is also a lower limit as 
we have not accounted for the loss of ions due to charge- 
exchange as they travel through the exosphere. This is, 
however, a negligible effect. Assuming a path through 
the magnetosheath in the GSE-XY plane from the sub¬ 
solar point around the Earth to roughly 120° from the 
sub-solar point (a path-length of ~ 55 Re) and the 0 +7 
charge-exchange cross-sec tion for a slow solar wind from 
Koutroumpa et al. (2006) (<r = 3.4 x 10 -15 cm 2 ), we find 


that less than a tenth of a percent of the ions recombine, 
and we have not even corrected this value for the number 
of 0 +8 ions that will have charge-exchanged to form 0 +7 
ions. Thus, correction of the ion density for the effects 
of charge-exchange is inconsequential. 


4.2. On the relative contributions of the magnetosheath 
and the heliosphere to SWCX “contamination” 

Efforts to understand the SWCX in the ChandraXMM- 
Newton era have been primarily driven by the problem of 
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Figure 12. A com parison of the magnetosheath size between the 
static|Spreiter et al.|(1966| model and the BATS-R-US model. The 
magnetopause distances are plotted with black crosses while the 
bowshock distances are plotted by red boxes. These were calculated 
from a subsample of the BATS-R-US simulations covering typical 
solar wind conditions. In all cases the distances are measured from 
the Earth along the GSE-X axis. 


removing the SWCX emission from obse rvations of dif¬ 
fuse emission. Kuntz & Snowden (2008) used multiple 
XMM-Newton observations with the same pointing but 
different epochs and observing geometries to show that 
SWCX emission is not solely a function of the local solar 
wind flux; some observations with low local solar wind 
flux showed SWCX emission while the majority of the 
observations at those solar wind fluxes did not. Further, 
the observations that showed SWCX emission when the 
local solar wind flux was low often had lines of sight 
that avoided the high emissivity regions of the magne¬ 
tosheath. However, they found no high local solar wind 
flux observations (nv > 10 9 cm -2 s -1 ) without SWCX 
emission. Finally, there were a number of observations 
that should have had lines of sight passing through the 
highest emissivity regions of the magnetosheath which 
showed no elevated SWC X emission, while other such 
observations did. Thus Kuntz & Snowden (2008) con¬ 
cluded that high local solar wind flux was necessary but 
not sufficient to predict SWCX emission. Although that 
SWCX emission observed when the local solar wind was 
low was probably due to local heliospheric emission, they 
c ould not exclude ma g netosp heric production. 


Carter & Sembay (2008) and Carter et al. (2011) 
combed the XMM-Newton archive for observations at- 
fected by SWCX. They searched for observations where 
the 0.5-0.7 keV band containing the oxygen lines showed 
variation not seen in a 2.5-5.0 keV “continuum ”_ban d. 
Their results were similar to those of IKuntz fe S nowden 
(2008j), finding that observations containing SWCX emis¬ 
sion tended to occur when XMM-Newton was between 
the Earth and the Sun, but were not limited to such 
lines of sight. For each observation they also compared 
the expected magnetosheath flux, as calculated from a 
static model of the magnetosheath and the ACE solar 
wind data. They found very poor correlation between 
the two values, suggesting both magnetosheath and he- 


liosphe ric co ntribu tions. 

Both Ku ntz fc S nowden] ([2008) and 
used the Spreiter et al. lJl9E 


Carter et al. (2011) 
H’or the size of the 


mode 

magnetoslieath given the solar wind density, speed, and 
temperature. We have compared the ma gnetopause dis 


tance and the bowshock distance of the Sp reiter et al. 
(1966) model with that found in the BATS-R-US simu¬ 
lations. The bowshock distance in the BATS-R-US sim¬ 
ulations was set to be the radius along the GSE-X axis at 
which the pressure drops to less than a tenth of the maxi¬ 
mum pressure. The magnetopause distance was set to be 
the maximum radius along the GSE-X axis at which the 
mag netic field li nes a re closed. Figure EH compares the 


Spreite r et al. (1966) and BATS-R-US values; the Spre- 


iter et al. ( |l966| ) consistently underestimates the mag¬ 


netopause and bowshock distances compared to those 
found by BATS-R-US . In the case of the bowshock, the 
underestimate is often 2-3 R#, while the magnetopause 
is usually off by ~ 1 Re - T hus, it is likely tha t both 


Kuntz & Snowden| (2008) and Carter et al. (2011) found 
a poor correlation between lines of sight passing through 
the nose of the magnetosheath and SWCX contamination 
because the lines of sight were actually passi ng several 
Re behind the regio n of substantial emission. Hen ley fc| 


Shelton (2010, 2012) extracted “blank sky” spectra from 
all possible XMM-Newton observations in order to study 
the Galactic halo. To study the SWCX emission they 
considered the oxygen line fluxes for cases of multiple 
observations of the same target. They looked, in partic¬ 
ular, at the difference between the flux of any particu¬ 
lar observation and the minimum flux for that direction. 
This quantity showed little correlation with the Q from 
the magnetosheath, _m part because their model of the 
magnetosheath ( S preiter et al.| 1966) had a stand-off dis¬ 
tance fixed to 10 Re without regard to the solar wind 
flux; their analysis suffere d from a more extreme f orm of 


the problem suffered by Kuntz & Snowden (2008 
Thus we conclude that many studies seeking SWCX 
emission from the sub-solar point of the magnetosheath 
have been looking in the wrong location. Whether the 
enhanced SWCX emission is observed when looking at 
the correct location has not yet been determined; a study 
of SWCX emission in the XMM-Newton archive using 
MHD models is currently ongoing at the University of 
Leicester. 


Henley & Shelton (2010 2012) also found “no univer¬ 
sal association between enhanced SWCX emission and 
increased solar wind flux”, a finding consistent with pre¬ 
vious studies. Given the lack of correlation between the 
local solar wind flux and Qh shown in §3.3, this is not 
surprising. (Note that XMM-Newton observes only a 
slightly larger region of the sky than was accessible to 
ROSAT .) However, since we have shown that the magne¬ 
tosheath emission produces a strong correlation between 
the | keV LTE flux and the local solar wind, but does 
not produce as strong a correlation between the | keV 


LTE flux (dominated by the emission from the O VII and 
O VIII lines) and the local solar wind flux, then the | 
keV LTE flux must be governed by some additional pa¬ 
rameter, most likely the oxygen ion abundances in the 
solar wind. If the ion abundance variation is enough to 
remove the bulk of the correlation seen from the mag¬ 
netosheath, then one would expect even less correlation 
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with local heliospheric emission. Unfortunately, the 0 +8 
abundance is poorly measured by ACE , and the 0 +7 is 
not much better measured, so it will be difficult to study 
this issue. 

In summary, the question of the relative contributions 
of the magnetosheath and the heliosphere, in general, to 
SWCX emission events is ill-posed, and studies in the 
oxygen lines are problematic. There should be a strong 
correlation between the local solar wind flux and mag¬ 
netosheath SWCX emission, even if previous attempts 
to verify that correlation failed due to incorrect expecta¬ 
tions for the location of the magnetosheath, and the use 
of the intrinsically more poorly correlated oxygen lines. 
Such correlation is not expected, in general, for helio¬ 
spheric SWCX emission. 


4.3. On SWCX emission prediction and SWCX 
emission removal 

Since the discov ery of the Hubble Deep Field SWCX 
emission event by Snowden et al. (2004), significant ef¬ 
fort has been placed on SWCX emission identification 
in | keV band observations, either to remove SWCX 
“contamination” for astrophysical reasons, or to iden¬ 
tify observations “b lessed” with SWCX e missio n for he¬ 
liospheric reasons. Kuntz & Snowden (2008) thought 
it a lost cause to determine whether any arbitrary ob¬ 
servation were effected by SWCX witho ut repeate d ob¬ 
servat ions of the same field, while [C arte r & Sembay 
(2008) showed if an observation was sufficiently long, 
then SWCX affected intervals c ould be identified through 
differential light-curve analysis. Henley & Shelton (2012) 
applied a local solar wind flux limit to remove SWCA 
emission affected observations. However, given the corre¬ 
lation analysis for the part of the sky observed by XMM- 
Newton , and excepting CME issues, this limit is as likely 
to exclude data without significant heliospheric SWCX 
emission enhancements as it is to include data with sig¬ 
nificant heliospheric SWCX emission enhancements. 

A second method to deal with SWCX emission is 
ex post facto modeling of the SWCX emission. Ef¬ 
forts to characterize the emission from the LHB have 
concentrated on repeated observations of nearby X-ray 
dark clouds to model and remove the SWCX emis¬ 


sion (e.g. Koutroum pa et al. 2011 Gupta et al. 2009). 
Koutroumpa ( 2012]) has shown, m particular, that mod¬ 
eling" the SWCX emission for a line of sight falling en¬ 
tirely within the equatorial flow yields reasonably reliable 
results, while modeling the SWCX emission for a line 
of sight dominated by the (poorly characterized) polar 
flow yields rather disappointing results. We note, how¬ 
ever, that these results were obtained at solar maximum 
when the distinction between polar and equatorial flows 
becomes problematic. At solar minimum, those flows 
are much more cleanly distinguished and the polar flow 
should be much more stable and amenable to modeling. 

A third method to deal with SWCX emission is to plan 
observations in order to minimize (or maximize!) its im¬ 
pact. Reference to Figure [To| suggests that to model the 
SWCX emission for an object within ~ 30° of the eclip¬ 
tic plane, one would want repeated observations of the 
object when it lies in the direction of the Parker spiral 
in order to optimize the correlation between the local 
solar wind flux and the SWCX emission that is actu¬ 
ally observed; observations in other directions will have 


significantly poorer correlations and thus poorer SWCX 
emission removal. Conversely, one might wish to ob¬ 
serve in such a way as to minimize the variation in the 
heliospheric flux. This would be particularly useful for 
shadowing observations where the shadow and the back¬ 
ground it is absorbing have to be observed separately. 
Figure [13] shows a map of the root mean square of the 
variation in Qh normalized by the mean of Qh • Al¬ 
though cfq h /Qh appears to vary more from simulation 
to simulation than does the correlation, some features 
are consistent. The highest variation is perpendicular to 
the Earth-Sun line in the same quadrant as the Parker 
spiral beyond the Earth. The opposite direction does 
not have an elevated variation. Thus it should be noted 
that high variability region is part of the sky accessible 
to XMM-Newton and Chandra , so SWCX variation can 
be reduced by specifying the part of the year in which 
the observations are to be made. 


4.4. On remote sensing of the magnetosheath 

One of the prime moti vations of this work is to pursue a 
suggestion first made in Collier et al. (2005), that SWCX 
emission from the magnetosheath could be used for re¬ 
mote sensing of the magnetosheath. This suggestion has 
spawned a number of mission prop osals to observe the 
magnetosheath, including STO RM (|Colhe r_et aL][2012) 
and AXIOM ( Branduardi-Raymont et al.||2012|), as well 
as smalle r proposals t o ob serve the cusps (Branduardi- 
Raymont & Mercier 2009). This is a particularly well 
posed problem as the SWCX production factor for hy¬ 
drogen can be derived directly from measurements of the 
object that we desire to study, albeit in the flanks of the 
magnetosheath rather than in the much brighter region 
of interest. 

As shown in Figure [l] the magnetosheath “breathes” in 
and out as the solar wind flux varies. Figure [14] demon¬ 
strates what would be seen with an X-ray imager 30 R# 
from the Earth, along the GSE-Y axis, looking back at 
the nose of the magnetosheath. The magnetopause and 
bowshock are clearly defined, as is emission beyond the 
bowshock due to the exosphere interacting with the free- 
flowing solar wind. There is a large literature exploring 
a wide variety of MHD phenomena that should occur 
within the magnetosheath, such as Kelvin-Helmholtz in¬ 
stabilities and flu x-transfer ev ents which may be observ¬ 
able in the X-ray (Su n et al.|2015 ). Whether these struc¬ 
tures can be observed depends critically on the count- 
rate and cadence of such an imager. 

As a first-order demonstration of the potential util¬ 
ity of such an instrument, we can determine the flux 
from the magnetosheath and compare it to that of the 
soft X-ray background. We have calculated the magne¬ 
tosheath emission as seen from a GSE coordinate (10,- 
30,0) R#, which produces images such as that seen in 
Figure [l4j For a characteristic value we calculated the 
emission along a line of sight through the peak emis¬ 
sion on the projected GSE-X axis. On the GSE-X axis 
the emission at the bowshock is roughly 20% of that of 
the GSE-X peak flux. We note that the magnetosheath 
closer to the dominant cusp tends to be even brighter 
than the “subsolar” flux. We have calculated the charac¬ 
teristic value for a number of simulations and have plot¬ 
ted the result in Figure [15) These simulation values can 
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Figure 13. The root mean square of the heliospheric SWCX emission meas ured over the 559 simulation steps divided by the mean of 
the heliosperic emission. The coordinate system is the same as in Figure [To| This map is the mean of four maps made from the same 
ENLIL simulation at 90° intervals around the Earth’s orbit. Although the individual maps are similar to the mean shown here, there is 
substantial variation. The direction of the He focussing cone has particularly large variation. That signal, which is at a different location 
in each individual map, was removed before forming this mean map. 



Figure 14. The 1 keV flux from the magnetosheath as viewed 
from GSE coordinate (0,-30,0). The Sun is to the right. The 
green circle to the left is the Earth. The white dashed contour 
is the projection of the bowshock. The white contours trace the 
emission in units of the typical background emission (2.2 count s —1 
deg -2 ). The straight white line is the projection of the GSE-X axis. 
The simulation is for early 2000-07-15 14:25:00, near the solstice, 
riswpVswp = 1-69 X 10 8 cm -2 s -1 , B^=-3.97 nT. The simulation 
was made with 0?5 X 0?5 pixels and is 41? 5 X 44°. No background 
has been added to this image. 

be compared the background. 

As Figure [l6| shows, the mean \ keV band (R12 band) 
flux over the entire sky is 6.16 x 10 -4 count s -1 arcmin -2 
(2.22 count s -1 degree -2 ) while the mode is 3.75 x 10 -4 
count s -1 arcmin -2 (1.35 count s -1 degree -2 ). These 
background values are significantly below the “subsolar” 
flux for all but the very calmest solar winds. Thus, imag¬ 
ing the magnetosheath is eminently feasible, given an in¬ 
strument with a large field of view and modest collecting 
area. 

5. SUMMARY 

In this work we have calculated the X-ray production 
factor for solar wind charge-exchange with H. The value 


Figure 15. The magnetosheath emission from a line of sight 
through the peak emission on the (projected) GSE-X axis as a 
function of the solar wind flux. The horizontal line is the mean | 
keV background level. 

is somewhat higher than theoretically calculated values, 
perhaps because the calculated values did not include all 
of the faint lines contributing to the emission in the | keV 
band. This value is rather robust, being based on the ra¬ 
tio of observed flux to the Q derived from BATS-R-US 
simulations; it is particularly robust for use with other 
BATS-R-US simulations being sampled for different ob¬ 
servation geometries. As a result, we have shown that 
remote sensing of the magnetosheath from high Earth 
orbit, looking back at the magnetosheath, is particularly 
promising. 

On the way to deriving the production factor, we 
may have resolved a few issues concerning SWCX emis¬ 
sion. We have shown how previous models of the mag- 
netosheath SWCX emission are unlikely to agree with 
XMM-Newton observations, simply due to the magne¬ 
tosheath not being in the expected location. We have 
also shown why analyses correlating SWCX emission in 
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thor and many of the coauthors were attending the Inter¬ 
national Space Science Institute (ISSI) workshop on solar 
wind charge-exchange in October 2013 in Bern, Switzer¬ 
land. Many of the coauthors also attended the IS SI work¬ 
shop on solar wind charge-exchange in January 2013. 

We are deeply grateful to IS SI for these very productive 
workshops. Our thanks also to the Community Coordi¬ 
nated Modeling Center at Goddard Space Flight Center 
for producing the extended ENLIL run, and for providing 
the large number of BATS-R-US runs used for this work. 
The BATS-R-US and ENLIL simulation results were pro¬ 
vided by the CCMC through their public “Runs on Re¬ 
quest” system (http://ccmc.gsfc.nasa.gov). The CCMC 
is a multi-agency partnership between NASA, AFMC, 
AFOSR, AFRL, AFWA, NOAA, NSF and ONR. 

Even with the large number of authors, a study of this 
type does not occur in a vacuum. We would like to thank 
all those who, through discussion, argument, and spirited 
altercations, helped make this a better paper. 
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